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Abstract. Low temperature and low frequency properties of a spin-boson model are investigated 
within a super operator and Liouville space formulation. The leading contributions are identified 
with the help of projection operators projecting onto the equilibrium state. The quantities of interest 
are expressed in terms of weighted bath propagators and static linear and nonlinear susceptibilities. 
In particular the generalized Shiba relation and Wilson ratio are recovered. 

PACS. 05.30. -d Quantum statistical mechanics - 66.35. +a Quantum tunneling of defects 



1 Introduction 

The spin-boson model is one of the simplest models to 
study the interplay of quantum mechanics and dissipation 
due to the interaction with an environment. It consists in 
a quantum mechanical spin j weakly coupled to a macro- 
scopic heat bath made up of harmonic oscillators. This 
and similar models of dissipative quantum systems have 
been widely studied in the literature [1,2,3]. 

The total system including the bath is described by a 
Hcrmitcan Hamiltonian and the corresponding Schroding- 
er or Heisenberg equation. The goal is, however, to obtain 

■ an effective description of the dynamics of the subsystem 
having averaged over the bath degrees of freedom. An ex- 
ample is the dynamics of a spin \ ruled by Bloch equa- 
tions. 

Especially for low frequencies and low temperatures 

■ one expects that an effective description might exist which 
is of the same structure as perturbation theory but with 
renormalizcd parameters. This is expected in analogy to 
the Landau-Fermi-liquid theory where effective masses 
and interactions can be obtained from static expectation 
values or experiments. 

For the spin-boson model results of this kind are 
known. For instance the ratio of the specific heat at low 
temperature and the static transverse susceptibility is de- 
termined by a bath propagator weighted with the spin- 
bath interaction. This has been found by Sassetti and 
Weiss [4] as generalization of the corresponding ratio for 
the Kondo problem, known as Wilson ratio. Perturbation 
theory yields the same result but with the transverse sus- 
ceptibility replaced by its free value. Another example is 
the Shiba relation for the Anderson model [5] general- 
ized again to the spin-boson model [4]. It states that the 



long time behavior of the transverse correlation function 
at zero temperature is again given by the static transverse 
susceptibility and the long time property of a weighted 
bath propagator. Again this result agrees with perturba- 
tion theory provided the susceptibility is replaced by its 
free value. 

The present paper deals with the low frequency and 
low temperature properties of a spin-boson system. In par- 
ticular a setup is investigated which allows to express the 
low frequency and low temperature properties in terms of 
bath propagators and static linear and nonlinear suscepti- 
bilities. This setup is based on a Liouville space formula- 
tion [6,7]. The fist section gives a brief survey of this for- 
malism. The spin-boson model is introduced in section 2 
and the weighted bath propagators are discussed. Section 
3 deals with the long time behavior of transverse and longi- 
tudinal correlation and response functions. This includes 
the Shiba relation and fluctuation dissipation theorems. 
The specific heat at low temperature and the Wilson ratio 
are discussed in section 4. As mentioned above the present 
analysis is based on the existence of certain static linear 
and nonlinear susceptibilities. The question under which 
conditions they do exist or do not exist [3], for instance 
because of localization, is not addressed. Furthermore the 
susceptibilities may depend on the high frequency cutoff of 
the bath oscillators resulting in non universal prefactors. 



2 Liouville space and super operators 

A convenient starting point is a formulation in Liouville 
space [6,7]. This is the linear space spanned by quantum 
mechanical Hermitean operators and each operator A is 
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considered as a vector \A) in this space. Super operators 
O are introduced as mappings of the quantum mechanical 
operators onto themselves. An example is the von Neu- 
mann equation for the statistical operator p(t) 



d^ (<) 



H(t),p(t) &C(t) p(t) 



(1) 



defining the Liouvillc super operator C{t). Here and in the 
following h = 1 is assumed. 

There are several possibilities to define a scalar prod- 
uct. In the following 



(A\B) = Tr AB 



(2) 



is used. 

The temporal evolution is described by the super op- 
erator U{t, t') which obeys 



-U(t,t') = C(t)U(t,t') 



(3) 



with initial condition U(t, t) = I and I\A) = \A). 

The measurement of an observable A can also be rep- 
resented as the action of a super operator A. with 



A\X) + 



(4) 



where the identity vector |1) obeys A\l) = \A). This al- 
lows to write the equilibrium correlation functions for time 
independent C and t > t' as 

C AB (t,t') = (l\AU(t,t')B\p). (5) 

The equilibrium statistical operator obeys C\p) = 0. 

Adding external fields to the Hamiltonian H(t) — ► H — 
hs(t) B response functions are defined as 



G AB (t,t') 



with 



5h B {t') ' 
B \X) & i 



(l\AU(t,t')B\p) (6) 



B,X 



(7) 



For spin | a complete set of quantum mechanical op- 
erators is formed by the identity and the Pauli matrices. 
They can be used as basis vectors of the corresponding 
Liouville space: 
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\y) ^ 
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\x) 
\z) 
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1 

-1 



(8) 



This results in a 4-dimensional representation 
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(9) 



Super operators are then 4x4 matrices, for instance 
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.(10) 



For particles, especially the bath oscillators, the Wig- 
ncr representation is used. The corresponding Liouvillc 
space is spanned by functions \X) X(p,x) of the coor- 
dinates x and momenta p. The scalar product is 



(Y\X) = J ApAxY{p 1 x)X{p 1 x). 



(11) 



The super operators for the measurement of position and 
momentum and the corresponding response operators are 



x <^=> x 



p <J=> p 



X <F> 



p ^ 



_d_ 
dp 
d_ 

dx 



(12) 



A general super operator O in Wigner representation is a 
function of four variables 0(x,p, x',p'). 

The Liouville space of the combined system and bath 
oscillators is the direct product of the individual spaces, 
for instance the four dimensional space of a spin \ and a 
Wigner representation for each bath oscillator. 



3 Spin-boson problem 

The Hamiltonian of a tunneling center or a spin i inter- 
acting with a thermal bath of N harmonic oscillators is 



H = -hAa x - 



■,ea 2 



Ttt T~] A k x k a z 

} 



(13) 



where A is the bare tunneling splitting and e a bias or 
anisotropy. Eventually the limit N — > oo is considered 
and the 1 / \/N scaling of the spin-bath coupling constants 
is explicitly written down. In the following the isotropic 
case e = is considered primarily, but generalizations to 
the general case are possible. 

The resulting Liouville operator (1) is 



C = \Ab x + \err z + -== A k \v z x k + x k & z } 



Pk P k +u 2 k x k x k 



} 



(14) 
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The bath is characterized by its density of states 



For the power law density of states (15) one finds 



J(w) = ^-V ^S(lj - Lu k ) = a0 1 - s w s e- a, / e (15) 

where the special case of a power law dependence is of in- 
terest. The case s = 1 is usually referred to as Ohmic case. 
A tunneling center in a crystal or glass requires a super 
Ohmic spectrum with s = 3. There is a high frequency 
cutoff representing the Debye frequency and the same 
frequency is used to define the dimensionless coupling con- 
stant a for s ^ 1. 

Because of the 1 / y/N dependence of the couplings the 
correlation and response functions of the bath oscillators 
are unperturbed. For an oscillator with frequency ujk the 
Xfe-Xfe-response function for t > is 



Fk(t) = — sin uo k t 



(16) 



and the corresponding correlation function at temperature 
T=l//3 is 



l~>k{t) = COSWfct 



with 



(17) 
(18) 



n(ui) — \ coth ^/3u). 

The action of the bath on the spin involves the weighted 
bath propagators 



and 



k 



Introducing the Fourier transforms 



-f 

Jo 



&te wt F{t) 



and 



/OO 
dte iwt D(t) 
-OO 

one finds for the imaginary part of F(uj) 



QF(u) = -si&i(u)J(\u\) 



(19) 



(20) 



(21) 



(22) 



(23) 



and for D(ui) 

D(lu) = im(\u}\)J(\w\) = 2n(u)3F(w). (24) 

This last expression is a special case of the quantum 
mechanical fluctuation dissipation theorem (FDT) [6,7] 



C AB (u) = 2n(u;)ZG AB {Lo) 



(25) 



which holds for any pair of operators A and B having the 
same parity with respect to time reversal. 



F(t) = a 2 r(s + 1) sin ((s + 1) arctan(<9i)) 



(l + (0i) 2 ) 



-(s+l)/2 



(26) 



and for T = 



D (t) = \a0 2 r(s + 1) cos ({s + 1) arctan(6>i)) 



x(l + (©i) : 



-(s+l)/2 



(27) 



The leading correction for finite temperature is SD(t) 
D(t) - D (t) with 



SD(t) = f 
Jo 



3 , J(P) 
auj— coswt 

= aT 2 (0/T) 1 - s d(Tt) 



(28) 



and 



d(r) 



COS XT 



/■OO 

= / dx — 
Jo 

OO 

= r ( s + X ) cos (( s + !) arctan(r/^)^) 



£=1 

x 



-(s+l)/2 



(^ 2 + r 2 ) 

r(s + l)C(s + 1) - |r(fl + 3)C(s + 3)r 2 
T(s) cos(s7r/2) t~ s + |r(s + 1) 



x sin(,s7r/2) r 



-s-1 



(29) 



A formal perturbation theory for expressions like (5) 
or (6) can be set up by writing 



U(t,t')=T 



exp 



(30) 



where T denotes appropriate time ordering. For the equi- 
librium statistical operator \p) = U(t',—oo)\i) is used. 
The initial state is not relevant assuming that the system 
has equilibrated at time t'. As a next step the exponential 
is expanded in a sum of products of free evolution oper- 
ators in spin space and bath operators. Finally all bath 
operators have to be contracted in pairs to bath response 
functions F(t) or correlation functions D{t). 



4 Long time behavior 

The following analysis of the long time behavior is based 
on two arguments: 

1) The system tends towards equilibrium 



U(t)-^\p){\\ 

t— »oo 



(31) 



2) The long time properties are ruled by the long time 
properties of the averaged bath response function F(t), 
(26), and the bath correlation function D(t), (27) and (28). 
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In order to isolate the leading contributions at long 
time a set of projection operators is introduced 



V = \P){1\ 

Ti = J2 {*fc|p)(l|P fc - P fc |p)(l|*fc} 



(32) 



XkPk - Pk X k 



}]|P)( 



II 



where the normal product : • • • : means that all bath re- 
sponse operators x k and p k have to be arranged to the 
left of |p)(l| and the operators x k and p k to the right. 
Furthermore for I > 1 each value of k may appear only 
once. 

Assume a projector Vi is inserted into the time evo- 
lution U(t,t') at some time t > t > t' 



U(t,t') ^U(t,t )V e U(t Q ,t'). 



(33) 



If the contraction of pairs of bath operators to bath propa- 
gators is performed, as outlined at the end of the previous 
section, exactly t such propagators reach from some time 
> to to some time < t . The proposal is that the lead- 
ing contribution for long time is given by the lowest order 
projector which can be inserted. 

Let me start with the investigation of the response 
function 

G zz (t) = (l\<r z U(t,0)& z \p) 

= {l\<r z U(t,0)& z U{0,-oc)\i). (34) 

Trying to insert Vo at some time t > to > results 
in vanishing contributions because (l\U(to,0) = (1| and 
(1\(T Z = 0. Next Vi is inserted. This yields 



(1| * z U(t,t )ViU(t ,0)& z U(0,-oo)\i) 

ft pt 
= I ds 
it Jo 

x(l\<r z U(t o ,0)v z \p) 

rt rto 



f ds [ ° ds'(l\<T z U(t,t )& z \p)F(s - s') 

Jtn JO 



= f ds [° ds'G zz (t-s)F(s-s')G zz (s'). (35) 

Jto J0 

It is proposed that (35) yields the leading contribution 
to G zz (t) for t — > oo and choosing for instance to = t/2 the 
s- and s'-integration can be performed replacing F(s — s') 
by F(t). This results in 



G zz (t) 



X 2 zz F(t) 

aQ 1 - s r{s + 1) cos(is7t)xL t' 1 - 



where 



Xz 



J de 



dtG zz (t) 



(36) 



(37) 



is static susceptibility which is assumed to be finite. This 
means that the asymptotic behavior of the time dependent 
transverse response function is determined by the static 
susceptibility and the averaged bath propagator. The bare 



parameters A, e and temperature enter only implicitly via 
X zz - For the Ohmic case s = 1 the first line of (36) is still 
valid but in the second line F(t) — > 2aO~ 1 t~ 3 has to be 
inserted. 

There are several points which have to be checked: In 
the integrals of (35) we have replaced F(s — s') by F(t) 
assuming that the main contribution comes from t-s«i 
and s' <ti t. This is actually the case for a bath spectrum 
with s > 0. If the asymptotic contribution (36) is inserted 
for G zz (t — s) or G zz (s') in (35) corrections ~ t _1_2s are 
obtained. 

Instead of 'Pi some other projector T 3 1 could have 
been inserted. For a symmetric system with e = 
the next nonvanishing contribution is obtained from 
£ = 3. A corresponding evaluation yields a contribu- 
tion \x 2 zzzz F{t)D 2 {t) where \ zzzz = 8d 3 (a z ) /de 3 is a 
static nonlinear susceptibility. This contribution vanishes 
_ f -3(i+ s ) for t < /? or - TH- 1 - 33 for t > (3. For e ^ 
corrections ~ t~ 2 ( 1+;s ) or ~ Tt~ x ~ 2s result from insertion 
of TV In any case (36) yields the leading contribution. 

Further contributions arise from insertion of 
Q = I~ J £2 e Ve- They are determined by the decay 
towards equilibrium. Within perturbation theory or 
NIBA [1,2,3] this decay is exponential and as a result 
the leading contribution (36) is indeed due to the first 
nonvanishing insertion of the projector Ve with lowest 
t = 1. The above analysis yields the restriction s > 0. In 
addition the existence of the static susceptibility (37) has 
been assumed which may impose further restrictions [10] 
on s or a. 

The Fourier transform 



G zz {uj) 

can be written as 

G zz (uj) = 

uj^O 



f 

Jo 



dte luJt G zz (t) 



Xzz 
Xzz 



- S zz (uj) 

■ XzzS zz (uj)Xz 



(38) 

(39) 
(40) 



with a self energy S zz (lu) vanishing for uj — > 0. With (21) 
and (36) 

£ zz (lo)^F(uj)-F(0)~lu s . (41) 

The analysis of the correlation function C zz (t) follows 
similar lines. The Fourier transform 

/oo 
dte iut C zz (t) (42) 
-oo 

may be written as 

C zz {uj) = G zz (uj)r zz (u)C zz (-uj) (43) 

defining the self energy F zz (t) . Again the leading contribu- 
tion at low frequencies or long time is due to an insertion 
of T\. Depending on whether a bath response function 
F or a bath correlation function D crosses the insertion 
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point at to several contributions are found. In analogy to 
(35) 



C zz (t)= [ ds f ds'G zz (t-s)F(s-s')C zz (s') (44) 

J to J— oo 

+ f ds [°ds'G zz (t-s)D(s-s')G zz (-s'). 

J to J — oo 

The second term is identified as due to the long time or 
low frequency part of F zz (t) and therefore 



r M (w)— ♦dm. 

uj— *0 



(45) 



The bath propagators fulfill an FDT (25) and therefore 
the same holds for the self energy S zz (uj) and r zz (ui). As a 
consequence the correlation- and response functions C zz 
and G zz also fulfill the FDT. This is not surprising. It 
shows the consistency of the above analysis. 

The leading time dependence for t — > oo is in analogy 
to (36) 

C zz (t) x 2 zz D(t) (46) 

t— »oo 

and with (27) C zz (t) ~ f" 8 " 1 for t -> oo and t < /?. In 
the opposite limit t — > oo and t » /3 an asymptotic decay 
- Tt- S is found from (28). 

The above result is a generalization of the Shiba rela- 
tion originally proposed for the Kondo problem [5] and 
generalized to the spin boson problem by Sassetti and 
Weiss [4] and others [8,9] for T = 0. Note that their defi- 
nition of x zz differs by a factor 4 from the present one. 

The analysis of the long time properties of the longi- 
tudinal response function 



G xx {t) = {l\a x U{t)a x \p) 



(47) 



follows similar lines. The lowest order projector which can 
be inserted is Vi resulting in 



G xx (t) 



ds dr ds' 

t a Jt JO 



X 



dr' (a x (t)a z (s)a z (r)) 
F( S -s')D(r-r')(a z (s')a x (0)a z (r')) 



+F(s - s')F(r - r') (a z (s')a x (0)a z (r')) \ (48) 

where (• • •) are correlation- response functions of higher 
order. Investigating the contribution of the first term, 
F(s — s')D(r — r') is replaced by F(t)D{t) resulting in 



G xx (t) ► x xzz F{t)D(t) 

t — >oo 



where 



Xx 



= 4 



d 2 (a x ) 



,d 2 (a z 



(49) 



(50) 



de 2 dAde 
is a static nonlinear susceptibility. The second term in 
(48) can be shown to be subdominant for t — > oo. As 
a result again the leading contribution at long time has 
been expressed in terms of nonlinear susceptibilities and 
weighted bath propagators. 



The longitudinal correlation function is obtained by 
applying fluctuation dissipation theorems. The resulting 
asymptotic expression is 



C xx {t) 



2 Xxzz -C(^) 



(51) 



decaying - t- 2 - 2s for t < f3 and - T 2 t~ 2s for t > (3. The 
same asymptotic behavior for T = has been obtained by 
Guinea [11] for s=l. The result obtained by Lang et al. 
[12] can not be compared since it refers to the longitudinal 
correlation function involving polaron dressed operators 
whereas (51) involves the bare operator a z . 



5 Specific heat 

In this section it is shown that thermal properties at low 
temperature can also be expressed in terms of static sus- 
ceptibilities and bath propagators. The free energy of the 
system at temperature T = 1//3 is 

F T = --|lnTre- /3H -lnTre- /3J?B | (52) 

where H is given by (13) and Hb is the Hamiltonian of 
the bath without coupling to the spin. 
Let me investigate the derivative 



dF T 
dA k 



1 



= ~JdtA k (a x {t){D k {t)a z {0) + F k (t)a x (0)} 
= -1 Jl\t A k [D k {t)G zz (t) + F k (t)C zz (t)}. (53) 

Because of the 1 / \/N dependence of the spin-bath cou- 
plings the functions G zz (t) and C zz (t) can be treated as 
independent on A k for N — > oo. This means that F is a 
linear functional of A 2 Dk(t) and A\Fk{t) respectively. 

The free energy Ft depends on temperature only via 
Dk{t). In order to evaluate 8 Ft = Ft — F the excess 



SD k (t) = D k {t-T) - D k {t-0) 
1 1 

= a 7 COS UJ k t 

0J k C 3 ^ - 1 



(54) 



is introduced. At low temperature the leading contribution 
is linear in 5D k (t) and 



8 Ft 



1 f c 



dtSD k (t)G zz (t) (55) 



dt8D(t)G zz (t) 



with SD(t) given by (28). For T -» the integral in (55) 
can be evaluated by using <5D(0) instead of 8D(t) because 
this function actually varies on a time scale f3 = 1/T — ► oo. 
This results in 



8 Ft 



T-»0 



hXzz8D(0) 



(56) 



= -\xz Z aT 2 {e/Tf- s r{s + i)as + i). 
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The specific heat is obtained from 

d 2 SF T 

C =- T ^ ( 57 ) 

resulting in 

C=^XzzaT'G 1 - s a(a + l)r(a + l)C(a + l). (58) 

which is a generalization of the Wilson ratio [4]. Note 
again the difference in the definition of the transverse 
static susceptibility. 
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